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Expressions are developed which permit calculation of the temperature-dependent thermal con- 
ductivity of a cylindrical specimen in which heat is generated internally, e.g., by radioactive decay. 
The information needed consists of the experimentally determined temperature distributions on the 
surfaces of the cylinder, the heat flow through a central circular area at one end of the cylinder, and 
the rate of internal heat generation (which in general may be position-dependent). Numerical coefficients 
arc tabulated for the case of uniform internal heat generation. The application of this calculation pro- 
cedure to published methods of thermal conductivity determination is shown and an example is given. 
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1. Introduction 

Subsequent to the publication of a paper |1|' in 
which a mathematical analysis of heat How in a right 
circular cylinder with arbitrary temperature boundary 
conditions was made, we were requested by D. L. 
McElroy of the Oak Ridge National Laboratory to con- 
sider the problem of determining the thermal conduc- 
tivity of a cylinder in which internal heat generation 
is also present. Since most methods of experimentally 
measuring thermal conductivity are riot readily applica- 
ble to self-heating specimens, and since there is a real 
need for determining the thermal properties of radio- 
active materials or of materials which are being irradi- 
ated, we felt the present analysis to be of sufficient 
interest to warrant publication even in the absence of 
any experimental data. 

2. Mathematical Development 

Consider a homogeneous, isotropic, opaque, solid 
right circular cylinder of radius, b, and thickness, /, 
with position-dependent heat generation, W(r), per 
unit volume per unit time, where r is the position vec- 
tor. In general, the thermal conductivity of the cylinder 
material may vary with temperature. The steady-state 
heat flow equation in cylindrical coordinates is 
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where temperature above an arbitrary datum is de- 
noted by the symbol, i\ the temperature-dependent 
thermal conductivity by /.(?), and the radial, longi- 
tudinal, and angular coordinates by r, z, and ft, respec- 
tively. Equation (1) can be reduced to a simpler form 
by introduction of a new variable 1 , //, defined by the 
relation 
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where &* is the value of k{v) at v = 0. Making the 
stitution (2), (1) reduces to 
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For Dirichlet boundary conditions (i.e., the tempera- 
ture distribution on the surface is specified), the solu- 
tion of (3) can be written as 



u = y-\- w, 



y=y(r, z, 0)~ 



w = w(i\ z, 6) z 
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of this paper. 



In (5), n is the outward drawn normal to the surface, 
dS' is an element of area on the surface, and the 
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integration is taken over the entire surface of the 
cylinder. In (6), dr is an element of volume and the 
integration is taken over the whole volume of the 
cylinder. The quantity G(r; r') is a time-independent 
Green's function which we choose to express in the 
two alternative forms [2, 3, 4, 5, 6]: 
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where J m is the ordinary Bessel function of the first 
kind and order ra, a mn is the nth positive root of 

Jm(a»ni) = 0, z< is the lesser of z or z, and z> is the 
greater of z or z; 



G(r; r') = 



1 



X X 



7T/P 



2 ^ cos m(d-d') 



m = — 3C H = 1 



sin (-y-j sin 



A77TZ 



/,„ (^) K. 



W77T; 



■K m i— I m i—j- 



Imimrrjl) 



(8) 



where I m and /JC»i are the modified Bessel functions 
of the first and second kind, respectively, and order 
m, r< is the lesser of r' or r, and r> is the greater of 
r or r. 

The physical meaning of the two components of 
(4) is easily seen. If we set W(r f ) = for all r', we see 
that w vanishes so that u = y. Thus y is the potential 
distribution due to the boundary conditions in the 
absence of internal heat generation. If we set u(r') = 
for all r' on the surface of the cylinder, we see that 
y vanishes so that u = w. Thus w is the potential 
distribution due to the internal heat generation with 
the boundaries held at zero potential. 

Equations (7) and (8) are completely equivalent; 
however, in some cases one may be preferred over the 
other, either because the integrations indicated in (5) 
or (6) are easier or because the infinite series resulting 
from these integrations are more rapidly convergent. 
If the surface integration of (5) or the volume integra- 
tion of (6) cannot be obtained analytically, y(r) or w(r) 
can be determined by numerical integration. 

Since we are considering the case of a homogeneous, 
isotropic material, W(v') would usually be expected 
to be axially symmetric for materials which self-heat 
by radio-active decay. In general W(r') would not be 
axially symmetric for cylinders which are heated by 
absorption of energy from an external source, e.g., 



irradiation by neutrons. For the specific case of 
axial symmetry, we assume the following boundary 
conditions: 



z = O^r^b u=f(r) 

z=l O^r^b u = g(r) 



(9) 



O^z^l r=b 



u=f(b)+[g(b)-f(b)\-+h(z\ 



where f(r) and g(r) are arbitrary potential distributions 
on the flat surfaces of the disk, and where, in order to 
insure continuity, we shall require that the function 
h(z) vanishes at the flat surfaces of the disk, so that 
h(0) = h(l) = 0. We can substitute these boundary 
conditions into (5) and obtain 
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where a n =a 0n and the coefficients A„, B„, and Ca are 
given by 
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In obtaining (10) from (5), we elected to use (7) to rep- 
resent G(r; r') on the flat surfaces of the cylinder and 
to use (8) on the convex surface. Equation (10) can also 
be written down without recourse to Green's functions; 
this was done in [1J. 

In the specific case of axial symmetry, we can inte- 
grate over 6 in (6) and obtain 

w(r, z) = 2tt \ dz' dr'r'W{r\ z')G(r, z; r' , z'), (14) 
Jo Jo 

where the Green's function now takes the two alterna- 
tive forms: 
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2.1. Uniform Heat Generation 

For the special case where W(r)=Wo, a constant, 
substituting (15) into (14) and performing the indicated 
integrations yields 
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eq (17) can be written in the form 
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The alternative expression for zc(r, 2), derived from 
(16), is 
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which can easily be derived, eq (10) can be written 
in the form 
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Equations (19) and (22) could have been derived 
without recourse to Green's functions for this specific 
case where W{y)=Wq since the leading terms in 
these two equations are particular solutions to Pois- 
son's equation for this case and are sufficiently simple 
to be written by inspection. In general, however, this 
could not be done. In principle, eq (6) provides us 
with a solution for w{r) regardless of the particular 
form of W(r) and hence is completely general. 

A practical case in which W(v) is essentially constant 
(and therefore (19) and (22) are applicable) is that of 
a sample which self-heats by radioactive decay and 
for which the mean free path of a fission fragment is 
much less than the dimensions of the sample. 

3. Calculation of Thermal Conductivity 

The total heat flow through a circle of radius a. where 
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From (4), this can be written as 
where 
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From (24), (25), and (26) we obtain 
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where the notation ra=l, 3, 5, under the summation 

sign indicates that the series is to be summed over This equation gives the thermal conductivity in terms 

odd values of // only. Making use of the identity of measurable or calculable quantities. Q" is the meas- 
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ured total heat flow through the circle r^a at z = (). 
The quantity (dyldz) z = o, which is the average tem- 
perature gradient which would exist over the circle 
/ ^ a at 2 = due solely to the boundary conditions 
but with no internal heat generation, can be computed 
from the observed temperature distribution on the 
surface of the cylinder using (5) and (25). The quantity 
(dwldz)z = {), which is the average temperature gradient 
which would exist over the circle r^a at z=0 due 
solely to the internal heat generation but with the 
boundaries maintained at zero potential, can be com- 
puted from the (assumed known) rate of heat genera- 
tion, W(r), in the cylinder using (6) and (26). Note that 
iv(r) is inversely proportional to A * so that the quantity 
k*{dwldz)z = o on the right-hand side of (27) is actually 
independent of A*. 

If we again consider the specific case of axial sym- 
metry, (10) can be substituted into (25) to obtain 
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where A„. B n , and C k are given by (11), (12), and (13). 
In the case of constant internal heat generation, 
W(r)= W , and (26) becomes 



Q' = TT(rlW {) VL 



(29) 



where fi = Cl(alb; I lb) has two forms corresponding to 
(19) and (22), respectively: 
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For large values of ///>, (30) approaches 
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Values of O are given in table 3 of [1|. There is an 
error in table 3 of [1 J in the value tabulated for a/b — 0.5, 
which should be 1.96630. 

Using a digital computer, numerical values for £l 
have been calculated for a range of values of a/b and 
lib and are given in table 1. Values were calculated 
using both (30) and (31) and the resultant two sets of 
values agreed to the number of significant figures 
given in table 1. Representative values are plotted 
in figure 1. 

As discussed in [1], the thermal conductivity, A*, 
corresponds to the reference temperature, jT*, from 
which the variable, v, is measured. In principle, T* 
may be arbitrarily selected and the mathematics will 
yield the corresponding value of A*. Thermal con- 
ductivity values may have to be computed by iteration 
as discussed in sections 4.2 and 5.2 of [lj. 

4. Example 

The above analysis is applicable to the several 
methods of determining thermal conductivity which 



tt 



ft 



a V 

b : b 



b" b 



46 2 A Jiia n alb) t . / / \ 

— 7 > — 7~, tanh a n —; 

al & aUtian) \ 2b) 

Ii(nnall) 



1 8^ 

2 7T :] a 



2 



nHoimrbll) 



In the special case, I lb -» 0, (30) reduces to 
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If lib -* 0, (31) also reduces to ft = l fa. 

In the special case, a/&-» 0, the two forms of H are 
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FIGURE 1. The coefficient 11 for the case of uniform heat generation. 

(«3ttJ The values of li corresponding to intermediate vales of alb lie between the two curves 

shown. 
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were discussed in [1J and also to other steady-state 
methods utilizing Dirichlel boundary conditions in 
which the temperatures at all points on the boundary 
are specified. In | l|. we discussed methods of thermal j n wn j c . n 
conductivity determination in which the investigators 
had assumed parabolic radial potential distributions 
on the flat surfaces of a cylinder and a linear longi- 
tudinal potential on the convex surface. This is a case 
having radial symmetry so that (10) can be used with 
the boundary conditions (9), in which /(r), g(r), and and 
h(z) are given by 
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and 



g{r)=Y x +Ei 



h(z) = 
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where Y and Y\ are the potentials (corresponding to y) 
at the centers of the faces of the disk, and E {) and £, 
are the potential differences between the edge and the 
center of the disk at z = () and z = L respectively. For 
these boundary conditions, (28) becomes 



are factors which depend only on the geometry of the 
cylinder and on the fraction of the disk from which the 
heat flow is measured. Numerical values of ^ and ^, 
are given in tables 1 and 2 of [1| for a range of values 
of a/b and l/b. 

In |1| we did not discuss the effect of internal heat 
generation. For the special case of parabolic radial 
potential distributions on the flat surface of the cylin- 
der, a linear longitudinal potential on the convex 
surface, and uniform heat generation throughout the 
cylinder, (27) becomes 



Table 1. The coefficient 11 as a function o/l/b and d/b for uniform heat generation 
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().() 


0.10 


0.20 


0.30 


0.40 


0.50 


0.60 


0.70 


0.80 


0.90 


1.00 
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0.0 


0.50000 


0.50000 


0.50000 


0.50000 


0.50000 


0.50000 


0.50000 


0.50000 


0.50000 


0.50000 


0.50000 


0.0 


.1 


.50000 


.50000 


.50000 


.50000 


.50000 


.50000 


.50000 


.50000 


. 1.9993 


,49872 


. 17328 


.1 


.2 


.50000 


.50000 


.50000 


.50000 


.49998 


.49995 


.49980 


.49924 


.49700 


.48786 


.44742 


2 


.3 


.49991 


.19989 


.49985 


. 19974 


.49949 


.49894 


.49766 


.49466 


.48745 


.46973 


.42244 


.3 


.4 


.49891 


.49883 


.49854 


.49796 


.49689 


.49495 


.49138 


.48477 


.47234 


.44841 


.39836 


.4 


0.5 


.49535 


.49511 


.49437 


.49295 


.49057 


.48670 


.48044 


.47030 


.45368 


.42584 


.37523 


0.5 


.6 


.48796 


.48754 


.48623 


.48384 


.48003 


.47423 


.46554 


.45257 


.43303 


.40301 


.35309 


.6 


.7 


.47657 


.47597 


.47413 


.47084 


.46577 


.45836 


.44779 


.43281 


.41148 


.38051 


.33201 


.7 


.8 


.46181 


.46107 


.45878 


.45478 


.44873 


.44015 


.42828 


.41205 


.38980 


.35872 


.31205 


.8 


.9 


.44463 


.44378 


.44118 


.43666 


.42995 


.42059 


.40793 


.39105 


.36851 


.33789 


.29327 


.9 


1.0 


.42599 


.42507 


.42226 


.41742 


.41030 


.40050 


.38747 


.37038 


.34800 


.31819 


.27570 


1.0 


1.2 


.38740 


.38643 


.38346 


.37840 


.37103 


.36108 


.34809 


.33144 


.31014 


.28250 


.24421 


1.2 


1.4 


.35038 


.34942 


.34651 


.34156 


.33441 


.32483 


.31247 


.29681 


.27706 


.25180 


.21737 


1.4 


1.6 


.31690 


.31598 


.31322 


.30853 


.30180 


.29282 


.28130 


.26683 


.24870 


.22571 


.19469 


1.6 


1.8 


.28749 


.28664 


.28406 


.27969 


.27343 


.26511 


.25448 


.24117 


.22459 


.20366 


. 1 7558 


1.8 


2.0 


.26202 


.26123 


.25884 


.25479 


.24900 


.24132 


.23153 


.21930 


.20411 


.18499 


.15944 


2.0 


2.2 


.24006 


.23932 


.23711 


.23337 


.22802 


.22093 


.21190 


.20064 


. 1 8668 


.16914 


.14575 


2.2 


2.4 


.22111 


.22043 


.21838 


.21491 


.20996 


.20339 


.19505 


.18465 


.17176 


.15559 


.13406 


2.4 


2.6 


.20471 


.20408 


.20217 


.19895 


.19435 


.18825 


.18050 


. 1 7086 


.15891 


.14394 


.12401 


2.6 


2.8 


.19044 


. 1 8985 


.18807 


.18506 


.18077 


.17509 


.16788 


.15889 


.14777 


.13384 


.11530 


2.8 


3.0 


.17794 


.17739 


.17573 


. 1 7292 


.16890 


.16359 


.15684 


.14844 


.13801 


.12502 


.10770 


3.0 


3.5 


.15272 


. 1 5224 


.15081 


.14840 


.1 1105 


.14038 


.13458 


.12737 


.11844 


.10726 


.09240 


3.5 


4.0 


.13368 


.13326 


.13201 


.12990 


.12687 


.12288 


.11780 


.11148 


.10367 


.09388 


.08088 


4.0 


4.5 


.11884 


.11847 


. 1 1 736 


.11548 


.11279 


.10923 


.10472 


.09910 


.09216 


.08346 


.07190 


4.5 


5.0 


.10696 


. 1 0663 


.10562 


.10393 


.10151 


.09831 


.09425 


.08920 


.08294 


.07511 


.06471 


5.0 




0.0 


0.10 


0.20 


0.30 


0.40 


0.50 


0.60 


0.70 


0.80 


0.90 


1.00 
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(42) 



As a specific example of the use of (42), we consider 
the method described by Hoch et al. [9], in which the 
heat flux from the flat surfaces of the disk is calculated 
from the temperature and the total hemispherical 
emittance of the specimen surface, using the Stefan- 
Boltzmann radiation law. That is, the heat flow, Q" , is 
given by 



Q"- 



2 eoT 4 



(43) 



where e is the total hemispherical emittance, o~ is the 
Stefan-Boltzmann constant, and T 4 is the average value 
of T 4 on the surface over the circle of radius a. Sub- 
stitution of (43) into (42) results in 



/[eoT 4 -/r () a] 

Yi-Yo+E&x-EoVo' 



(44) 



With the simplified boundary conditions assumed by 
Hoch et al., namely Fi = Fo and E\ = Eo, this reduces 
to 



/[eoTwr ni 

E [Vi-Vo\ 



(45) 



If this is evaluated at the center of the disk (i.e., a = 0), 
as was done by Hoch et al., the thermal conductivity 
is given by 



** = 



l[ea% 4 -lW () £l\ 



(46) 



where T {) is the (absolute) temperature at the center of 
the disk at 2 = 0, and ^o, ^i and fl are to be evaluated 
at a/6 — 0. If there is no heat generation (i.e., Wq = 0), 
equation (46) reduces to a form which is exactly equiv- 
alent to eq (13) of Hoch et al. [9J. 

We postulate the following data as having been 
acquired experimentally: 



6=1.00 cm / = 1.80 cm 

e = 0.17 r () = 14.3 W/cm : * 

r = 2360°K Eo = 28 deg. 

From tables 1 and 2 of [1J for the case a/6 = 0.0 and 
//6 = 1.8, we obtain 

*, = 0.87364 ^o =~ 2.85220. 

From table 1 of the present paper, we obtain 

a = 0.28749. 

Using o~ = 5.6697X 10~ 12 W/cm- deg 4 , substitution of 
the above values into (46) yields 

k* = 0.388 W/cm deg 

as the thermal conductivity of the sample at the tem- 
perature corresponding to v=0. 
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